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FIELD OF THE INVENTION 

This invention relates to pulse oximeters which noninvasively measure the oxygen 
10 saturation of hemoglobin in arterial blood and, in particular, to an improved system for 
performing these calculations and for discriminating valid signals fi-om invalid data such as 
caused by probe-off conditions or motion artifacts. 

PROBLEM 

It is a problem in the field of pulse oximeters to accurately measure the oxygen 
15 saturation of the hemoglobin in arterial blood without having a significant error content due to 
ambient noise. It is also a problem to determine when the probe used to perform the 
measiu-ement is producing data with a noise component that is too large to provide an oxygen 
saturation reading of sufficient accuracy. The oxygen saturation (Sp02) of the hemoglobin m 
arterial blood is determined by the relative proportions of oxygenated hemoglobin and reduced 
20 hemoglobin in the arterial blood. A pulse oximeter calculates the Sp02 by measuring the 
difference in the absorption spectra of these two forms of hemoglobin. Reduced hemoglobin 
absorbs more light in the red band than does oxyhemoglobin while oxyhemoglobin absorbs 
more light in the inft-ared band than does reduced hemoglobin. The pulse oximeter includes a 
probe that is placed in contact with the skin, in a vascularized tissue region, either on a flat 
25 surface in the case of reflectance probes or across some appendage in the transmission case. 
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The probe contains two light emitting diodes (LEDs), each of which emits light at a specific 
wavelength, one in the red band and one in the infrared band. The amount of light received 
after interacting with the tissue adjacent to the LEDs is measured many times a second at both 
wavelengths. 

5 The tissue contains arterial, capillary and venous blood as well as muscle, connective 

tissue and bone. Therefore the red and infrared signals received fix)m the probe contain a non- 
pulsatile component which is influenced by the absorbency of tissue, venous blood, capillary 
blood, non-pulsatile arterial blood, the intensity of the light source and the sensitivity of the 
detector. The pulsatile component of the received signals is an indication of the expansion of 

10 the arteriolar bed with arterial blood. The amplitude of the pulsatile component is a very small 
percentage of the total signal amplitude and depends on the blood volume change per pulse and 
the Sp02. The received red and infrared signals have an exponential relationship to their 
respective incident intensities. Therefore, the argument of the received red and infrared signals 
have a linear relationship and these received signals can be filtered and mathematically 

15 processed using either derivatives or logarithms. The effects of different tissue thicknesses and 
skin pigmentation can be removed from the received signals by normalizing the processed 
signal by a term that is proportional to the non-pulsatile portion of the received signal intensity. 
Taking the ratio of the mathematically processed and normalized red and infrared signals results 
in a number which is theoretically a function of only the concentration of oxyhemoglobin and 

20 reduced hemoglobin in the arterial blood, provided that the concentration of dyshemoglobins in 
the arterial blood is sufficiently small. 

This data is significantly impacted by the presence of noise, which is manifested in many 
fohns. Aixy phenomena, vy^hether mechanical or electrical or optical, that causes an artifact in 
the pulsatile component of the received signals compromises instrument perfomiance. An 
25 example is that any transient change in the distance between the light emitting diodes and the 
detector can result in an error signal at both wavelengths. These pulse sources can course 
annoying false positive alarms which are of concern in a critical care setting because 
instruments with frequent alarms are often ignored or the alarms are disabled. Motion artifacts 
can be caused by patient movement and firequently mimic vascular beats with frequencies well 
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within normal physiological ranges. A second source of noise is the introduction of ambient 
light into the probe. Any light incident on the detector and not originating from the light 
emitting diodes is considered noise. Many of these noise sources are not easily filtered out of 
the resultant signal and represent a significant error component in existing pulse oximeters. 

5 SOLUTION 

The above-described problems are solved and a technical advance achieved in the field 
by the improved system for noninvasively calculating the oxygenation of hemoglobin in arterial 
blood using a pulse oximeter. This improved system takes advantage of the basic statistical 
property of pulse oximetry signals that the measured blood oxygen saturation appears as a 
10 constant over a small set of measurements. Properly processed sets of red and infixed signals 
have a linear relationship therebetween. 

The processing has several steps. First, the received red and infixed signals are 
collected from the probe detector and pre-processed. This pre-processing includes removal of 
ambient tight and prefiltering to remove noise. The pre-filtering can be a combination of linear 

1 S filtering to remove unwanted additive noise and non-linear filtering to remove noise spikes. The 
filtered red and infrared signals consist of both a small magnitude pulsatile component which 
carries information about the oxygen saturation of the hemoglobin, and a plurality of large 
magnitude non-pulsatile components. The filtered signals are mathematically processed using 
either derivatives or logarithms. The results of the mathematical processing are a set of red 

20 values which are directly proportional to the red optical Extinction, and a set of infrared values 
which are directly proportional to the infrared optical Extinction. The red and infixed signals 
travel through the same pulsatile path when the data is good, which leads to a technique for both 
identifying good data and extracting the ratio of red optical absorption to the infixed optical 
absorption for that set of data measurements. 

25 The oxygen saturation is essentially constant for a set of measurements taken over a 

short interval of time; e.g. less than a second. This implies that, when the data is good and the 
optical paths are the same for the red and infrared signals, the ratio of the mathematically 
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processed red signal to the mathematically processed infrared signal is also a constant, except 
for residual noise. Consequently, for good data, a plot of the simplified infrared data points 
versus the simplified red data points yields points that are tightly scattered around a straight line. 
A "best-fit" straight line can be computed for these data points using standard mathematical 
5 techniques such as linear regression. The slope of that best-fit-line is proportional to the RRatio, 
which is defined as the ratio of red optical absorption to the infrared optical absorption for that 
data set. The RRatio carries the desired information about oxygen saturation of hemoglobin in 
the arterial blood. 

The processed signals also contain additional inforaiation. If the plotted points are 
10 widely scattered about the best-fit-line, the excessive scatter indicates that there is excessive 
noise or probe artifacts in the received signals. This can be an indication the probe has fallen off 
the patient and/or there is patient motion. The degree of scatter of the points can be measured 
using standard mathematical techniques, such as the linear correlation coefficient. The measure 
of the scatter of the plotted points provides a quality measure of the data sample. Other 
1 5 statistical techniques could also be used to check the linearity of the data; for example, higher- 
order moments or higher-order fits. Still another test of data quality can be the intercept of the 
best-fit-line. The intercept is very close to zero for good data. The computed error measures 
can be compared to a failure threshold, possibly over a series of measurements. If there is too 
much error, then an alarm can be generated for the user. 

20 The minimum size of the measured data set required for high confidence Sp02 

calculations is small and therefore permits faster response of the pulse oximeter to changing 
values of the Sp02. It is not necessary to wait a pulse period to acquire enough samples, and 
high confidence values can be computed in less than a second. 
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BRIEF DESCRIPTION OF THE DRAWING 

Figure 1 illustrates in block diagram form the overall architecture of the improved 
system for Sp02 calculation; 

Figure 2 illustrates in block diagram form the functional elements used by this system to 
5 perform the calculations; 

Figures 3-7 illustrate additional details of the various processes used in the improved 
system for Sp02 calculation; 

Figure 8 illustrates in graphical form the various components of the input signal from the 
probe; and 

10 Figure 9 illustrates in graphical form an illustrative set of data which is used in this 

system. 

DETAILED DESCRIPTION 

This pulse oximeter system takes advantage of the basic statistical property of pulse 
oximetry signals that the measured blood oxygen saturation appears as a constant over a small 
15 set of measurements. Properly processed sets of red and infrared signals should have a linear 
relationship therebetween. Preprocessing of the sets of red and infrared signals received from 
the probe removes ambient light components and linear and non-linear prefiltering also removes 
noise to provide more accurate measurement data. 

The fihered signals are mathematically processed using either derivatives or logarithms. 
20 The results of the mathematical processing are a set of red values which are directly 
proportional to the red optical absorption, and a set of infrared values which are directly 
proportional to the infrared optical absorption. The presence of good data is indicated by a 
linear relationship between the mathematically processed red values and the mathematically 
processed infrared values. 
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There are several types of appropriate mathematical processing. For example, one 
method is based on a derivative technique which finds the numerical difference between two 
successive red measurements. This difference results in a red term that is proportional to the 
product of the red optical absorption and the pulsatile path, provided that the sample rate is rapid 
5 enough for the time between measurements to be very small. Taking a similar difference 
between two successive infrared measurements provides an infrared term that is proportional to 
the product of the infrared optical absorption and the pulsatile path. For good data, the pulsatile 
path components are identical for both red and infrared signals, so that the ratio of the 
differentiated red and infrared terms cancels out the path length. The effects of different tissue 
10 thicknesses and skin pigmentation are removed by normalizing the processed signal by 
removing the terms that are proportional to the red and infrared non-pulsatile incident 
intensities. 

An alternate method of mathematical processing is to use logarithms instead of 
derivatives. In this approach, flie logarithms of the ratio of two successive red measurements 
15 are taken, which is equivalent to finding the difference between the logarithms of two 
successive red measurements. The same process is appHed to the infrared signals. Then the 
ratio of the logarithmic differences of the red and infrared signals yields a term that is a ftinction 
of the red and infrared absorptions and is independent of path length. This technique directly 
compensates for the variations in thickness of different tissues and skin pigmentation. 

20 Regardless of whether the derivatives or logarithms are used for mathematical 

processing, the linearity of the mathematically processed red and infrared signals allows the use 
of statistical techniques such as linear regression and linear correlation to fit a straight line to the 
set of pairs of processed red and infrared data points and to measure the goodness of that fit. 
One result of this analysis is a linear regression slope and intercept and the other result is a 

25 goodness-of-fit correlation coefficient. The linear regression slope is a measure of the best fit 
straight line to the ratio of the red absorption coefficient and infrared absorption coefficient. 
This ratio is called a best-fit RRatio. This best-fit RRatio is then used in an empirical 
calibration formula to compute the best-fit Sp02 value. The minimiun size of the data set 
required for high confidence Sp02 calculations is relatively small, and permits faster response 
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of the pulse oximeter system to changing values of Sp02. It is not necessary to have a full pulse 
period, and responses under one second are feasible. 

The goodness-of-fit correlation coefficient is a measure of linearity between the 
mathematically processed red and infrared signals. If the correlation coefficient is. low, then the 
5 data is a poor fit. This can indicate that there is a probe-ofF condition or patient motion. By 
tracking a short history of computed correlation coefficients against pre-determined failure 
thresholds, error conditions can be detected and appropriate alarms generated for the user. A 
large magnitude of the computed intercept of the best-fit-line also indicates bad data. 

System Architecture 

10 Figure 1 illustrates in block diagram form the overall architecture of a typical pulse 

oximeter system including the apparatus of the present invention. The pulse oximeter system 
100 consists of a probe 101 connected to probe interface circuit 102 by means of a set of 
electrical conductors 103. The probe 101 consists of an exterior housing 104 that applies the 
active elements of the probe to the tissue under test, such as a finger 105, containing an arterial 

15 blood flow that is to be monitored. Included within housing 104 is a pair of light emitting diodes 
111, 1 12 and at least one corresponding light detector 113. The light emitting diodes 111, 112 
each produce an output light beam of predetermined wavelength which is directed at the tissue 
under test 105 enclosed by housing 104, The light detector 1 13 monitors the level of light that 
is transmitted through or reflected from the vascularized tissue 105. In order to distinguish 

20 between the light beam produced by the red 1 1 1 and infrared 112 light emitting diodes, these 
light emitting diodes 111, 112 are synchronously sampled. Ambient light is measured in 
housing 104 by light detector 113 making measurements when the light emitting diodes 111, 
112areoff. 

The signals received by probe interface circuit 102 from light, .detector 113 are analog 
25 signals and are typically processed by additional analog circuitry and then converted by an 
analog-to-digital converter circuit into sets of digital measurements which are stored in memory 
106. Data processing circuit 107 mathematically processes the digitized measurements stored in 
memory 106 to compute the oxygenation level of the hemoglobin in the arterial blood in tissue 
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105. It is also possible that analog circuitry can be used to perform some of the mathematical 
operations described herein as performed by data processing circuit 107, such as taking 
derivatives or logarithms. 

The data processing circuit 107 also computes ibt pulse rate, and a pulsitility index. The 
5 results of the oxygen saturation computation and pulse rate are displayed numerically (115) via 
display driver 109 and the associated display 1 14, 1 15 while the plethysmographic waveform is 
typically displayed graphically via display 114. The pulsitility index can be done with 
numerical or graphical methods. In addition, data processing circuit 107 detects anomalies in 
the input data stored in memory 1 06 and "sets error codes which may be used to" activate an 
10 alarm generation circuit 108 to produce an audible and/or visual alarm to a user to indicate a 
data error condition. The operation of data processing circuit 107 is disclosed in additional detail 
below and, for the purpose of this disclosure, it is assumed that the other elements disclosed in 
Figure 1 are the conventional components found in existing pulse oximeter systems. 

Signal Components 

15 Figure 8 illustrates in graphical form (not to scale) the various components of the signal 

produced by the light detector 1 13 as a result of a light beam interacting with vascularized tissue 
105. The light detector output signal consists of a large magnitude non^pulsatile component and 
a small magnitude pulsatile component. The non-pulsatile component represents light 
remaining after absorption due to a combination of venous blood flow, tissue, bone, and 

20 constant arterial blood flow while the small pulsatile component is caused by the light 
absorption due to pulsatile arterial blood flow that is to be measured. Since the light emitting 
diodes 1 1 1, 1 12 are sampled in rapid succession, the data signals produced by the light detector 
113 and transmitted to the probe interface 102 consist of a series of data points which are 
digitized and stored in memory 106. Since the red 1 11 and infrared 112 light emitting diodes 

25 are sampled in rapid succession, the digitized data points produced consist of a plurality of sets 
of measurements, with one set corresponding to samples of the red intensity, the other set 
corresponding to samples of the infrared intensity, and a third set corresponding to ambient 
light. 
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Ideally, the ratio of the normalized derivative (or logarithm) of the red intensity to the 
normalized derivative (or logarithm) of the infrared intensity is a constant. This constant is 
indicative of the partial oxygenation of the hemoglobin in the arterial blood flow. It is obvious 
that this ratio changes as Sp02 changes but, for a short interval with rapid enough sampling 
5 rate, the ratio remains constant. 

As noted above, the actual data received by the probe interface circuit 102 can include a 
fairly significant noise component which is caused by a number of sources including motion of 
the finger 105, the introduction of ambient light into the housing 101, and various sources of 
electrical noise. These noise components skew the values of either or both of the magnitudes 

10 measured in each set of data points destroying the correct relationship between the red and 
infirared signals. Existing pulse oximeter circuits make use of various filtering techniques to 
minimize the impact of the noise on the Sp02 value measured by the system. However, a 
severe limitation of existing filtering circuits is that very little indication is provided to the user 
of the severity of the noise introduced into the measurements. The accuracy of the readings 

15 produced by a typical pulse oximeter system can be severely compromised by the probe 101 
being improperly affixed to the patient's finger 105 without the user being able to detect this 
condition. In addition, only a single pass filtering technique is typically used to remove only the 
most egregious instances of noise in the system. The apparatus of the present invention 
overcomes the limitations of existing pulse oximeter systems by processing the received data in 

20 a manner that detects excessive noise introduction in the system and provides improved noise 
rejection in determining the Sp02 measured by probe 101. 

Mathematical Background 

In this system, the key parameter is called the RRatio. This parameter measures the ratio 
of the normalized derivative (or logarithm) of red intensity to the normalized derivative (or 
25 logarithm) of infrared intensity. For good data, this parameter corresponds to the ratio of the red 
arterial optical absorption to the infrared arterial optical absorption. The significance is the 
. RRatio can be understood by examining the optical behavior of light as it passes through tissue. 
The light is scattered and absorbed by all the tissues, but the light passing through a pulsing 
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artery or arterial bed will see a moving path length. The other tissues are unmoving and 
contribute to the steady non-pulsatile signal, but not to the time-varying pulsatile signals. The 
absorption of light by arterial blood is assumed to be only a function of the oxygenation state of 
the hemoglobin. Other basic assumptions are that the red and infrared light travels along 
5 essentially the same optical path, and that all the hardware circuits do not introduce any bias into 
the signal extraction. 

The waveform of the optical signal is determined by changes in the pulse length, Lp, and 
the basic optical signal is described by an exponential relationship according to Beer's law. 
This means that the intensity of red, iRcd, and the intensity of the InfraRed, Iir, light absorbed by 
10 the blood in a pulsing capillary bed can be described by the following equations: 

iRed = Ambient + Inpred * Exp[-ERed*C* Lp] (1) 

Iir = Ambient + Inpir * Exp[-E|K*C*Lp] (2) 

where the Red or IR optical absorption is given by the product of the Red Extinction 
Coefficient, ERed, or InfraRed Extinction Coefficient, Eir, and the hemoglobm concentrations, 
IS C. The concentration, C, and pulsatile path length, Lp, are the same for the red and infrared 
light. Inpred is the non-pulsatile red intensity of the signal and is proportional to the product of 
the light emitting diodes intensity with attenuation factors from propagation through non- 
pulsatile media and sensitivity and gain factors from the detection electronics. Inpir is the 
corresponding non-pulsatile infrared intensity. 

20 The pulsatile exponential function puts an amplitude modulation on the npn-pulsatile 

intensity. The equation for the non-pulsatile intensity Inp, has several terms. 

Inp = Io*A*Exp[-£(Ex*Cx*Lx)] (3) 

The term lo is the incident intensity. A is a composite scale factor caused by scattering losses, 
losses in detection, and gain in detection circuits, there is a composite exponential absorption 
25 as the light propagates through non-pulsatile intervening tissue, which can be described by a 
summation of the absorptions by each intervening tissue zone. 
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The objective of the following mathematical manipulation is to obtain the ratio of the 
Extinction coefficients of the Red, ERcd, and the InfraRed, Ejr, signals. This ratio is used with a 
calibration equation to determine how well the blood is oxygenated. It is assumed that the 
ambient light is perfectly subtracted from the received signals. 
5 It is necessary to compute the Extinction coefficients from the exponentials. This can be 

done with the use of either logarithms or derivatives. The preferred method disclosed herein 
uses derivatives. If the samples are closely spaced relative to the changes in pulsatile path 
length, Lp, then a series expansion of the exponential function can be used. This means that the 
exponential really looks like a constant plus a term that is linearly proportional to the path 
10 length. In other words^ only the first two terms of the binomial expansion need to be retained. 
Assume that the difference between two closely spaced in time non-pulsatile intensities is zero. 
Taking samples at time t2 and time tl, the derivative of the red mtensity can be approximated as 
the difference between the two samples: 

DRed = lNPRED*(Exp[-ERed*C*Lp(t2) ] - ExphERed*C*Lp(tl)])= 

15 Inpred* ([l-ERed*C*Lp(t2) ] - [l-ERea*C*Lp(ti) ] ) = 

Inpred *(ERed*C*[Lp(t2) - Lp(ti) ]) (4) 

In equation 4 Lp(t2) is the pulsatile path length at time 2 and Lp(tl) is the pulsatile path length at 
timel. 

The derivative of the infirared intensity can be approximated in a similar manner. 
20 DIR = Inpir* (Exp[-EiR*C*Lp(t2) ]-Exp[-E,R*C%(tl) ] ) - 

Inpir* ( [ l-E,R*C*Lp (t2) ] - [ l-EiR*C*Lp (tl) ] ) = 

. . Inpir* (EiR*C*[Lp(t2)-Lp(tl)]) (5) 

Note that both DRed and DIR are directly proportional to the Extinction coefficients as 
long as Inpred and Inpir are constant over that sampling interval. Note also that the term 
25 C [Lp(t2)-'Lp(tl)] is common to both DRed and DIR. This means that the ratio of DRed and 
DIR is independent of the concentration and path length changes. 
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DRed/DIR = (lNPRED/lNPiR)*(ERed/E,R) (6) 

NOW the factor which is due to the ratio of received non-pulsatile intensities can be 
removed in several ways. In one method, analog hardware is used to separate the received 
intensities, by finding the low-pass DC values of the Red and Infrared. 
5 (Inprei/Inpir) ~ (Inpred dc/Inpirdc) (7) 

In the preferred method, the sum of the two or more time samples is used 

Inpred = lNPRED(t2) + iNPRECKti ) (8) 

InPIR = lNPIR(t2) = lNPIR(tl ) (9) 

An instantaneous RRatio is found by multiplying the ratio of derivatives to the ratio of the 
10 incident intensities. 

instantaneous RRatio = DRed/DIR*(lNPiR/lNPRED) = (ERcd/EjR) (10) 
Key Equations 

The instantaneous RRatio is a number computed from only two measurement pairs of 
red and infrared values, which can be corrupted by noise. There is another way to computed an 

15 RRatio that forms a best-fit to a larger number of measurements pairs and reduces the sensitivity 
to noise. Note once again that the Extinction Coefficients and incident intensities are nearly 
constant provided that the sample rate is fast compared to any physiological changes. This 
means that valid data must obey the equation of a straight line. Consequently for a small set of 
closely spaced samples, it must be true that: 

20 DRed*lNPiR = SlopeA*(DIR*lNPRED) + noise, and (11) 

where 

SlopeA = ERed/EiR (12) 
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An alternate way to express the linear relationship may be more applicable in some 

cases. 

DRed - SlopeB*DIR+noise, and (13) 
where 

5 SlopeB=(lNPRED/lNPIR)*(ERed/EiR) (14) 

Equations 11 and 12 state that each value of DIR should be scaled by multiplying by 
Inpir and that each value of DRed should be scaled by multiplying by Inpred- Then a linear 
regression of all the scaled points in the set yields a line whose slope is the best fit to the ratio of 
the Red Extinction to the Infrared Extinction which is the desired RRatio. 
10 best-fit RRatio = SlopeA (15) 

If the incident intensities are constant over the data set, then it is not necessary to scale 
each value before doing the linear regression. That case corresponds to Equations 13 and 14, 
which indicate that the linear regression slope can be computed on DIR and DRed and then 
normalized by the ratio of the incident intensities. 
1 5 best-fit RRatio = SIopcB/CInpred/Inpir) (16) 

After finding the best-fit RRatio, the Sp02 is computed &om an empirical calibration formula. 

Equations U to 16 are the key equations. They state that there is a linear relationship 
between the DIR and DRed provided that the data is good. Consequently, the slope of the linear 
regression fit to a set of data samples yields a best-fit estimate for the ratio of the Extinction 
20 Coefficients, RRatio. The goodness of the fit can then be used as a test for checking the quality 
of the data. Any deviation fix)m the linearity relationship indicates a reduced quality of the data. 

For completeness, we show how the linear regression method can be used when the red 
and infi-ared signals are processed with logarithms instead of derivatives. In this case, define the 
difference between the logarithms of two subsequent red and infirared measurements. 
25 DLogRed = 

(log[lNPRcd{t2) ]-L0g[lNPRED(tl)l ) + ERed*C (Lp(arLp(ti)) (17) 
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and 

DLoglR = 

{log[lNPiR{t2)]-Log[lNPiR(tl)l) + E,R*C*(Lp(t2)]-Lp(ti)) (18) 

If two samples are closely spaced in time, the difference between two non-pulsatile intensities 
5 goes to zero. Then: 

DLogRed/DLogIR = ERcd/EiR (19) 
Since the Extinctions are constant, this means that: 

DLoglR = SlopeA*DLogRed + noise (20) 
where 

10 SlopeA = ERcd/EiR (12) 

- Equation 20 and Equation 1 1 have the same linear form. The resuUs from equation 20 and 1 1 
differ only if there are significant variations in incident intensity or large amounts of additive 
noise, or the sample spacing is too large for the series expansion to be vaUd. 

System Impiementation 

15 Figure 2 illustrates in block diagram form the various components used within the data 

processing circuit 107 of the pulse oximetry system of the present invention in order to provide 
improved data filtering and data acquisition. The analog data signals received from probe 101 
are filtered by analog hardware and digitized by an analog to digital converter in probe interface 
102 prior to application to preprocess input data from probe circuit 20 L The preprocessing 

20 circuit 201 includes prefilter circuit 301 (Figure 3) which uses digital filtering techniques for 
noise removal, fi^uency selection and data rate control. Nonlinear digital filters can be used in 
the prefilter circuit 301 of preprocessing circuit 201 to reduce ttie influence of input data that are 
statistically beyond the range of the other data values acquired by probe 101. These nonlinear 
digital filters can be implemented for example by median filters, trimmed mean filters and 

25 morphological filters. 
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The input data from probe 101 is typically decomposed into its non-pulsatile and 
pulsatile sub-elements in probe interface 102 in order to provide accurate measurements of these 
components. The pulsatile component typically represents 0.5% of the total input signal and the 
decomposition of the input signal into pulsatile and non-pulsatile components permits accurate 
5 analog to digital conversion of each of these components. Therefore, once these two sets of 
components have been filtered by the prefilter circuit 301 of preprocessing circuit 201, they are 
applied to reconstruct red, IR, and ambient circuit 302 where the components are used to return 
the acquired input signal back to its original form. The reconstructed infrared and red data 
signals are in digital form and the signal reconstructed within reconstruction circuit 302 is used 
10 for pulse computations and display 115. The red and/or infrared signals are further processed 
for waveform display 1 1 4 on the pulse oximeter system. 

For some hardware configurations, the reconstructed red and infi:ared signals are 
transmitted to process ambient input data circuit 202 which is comprised of the elements 
illustrated in Figure 4. The ambient signal is removed from the red and infrared signals in 

15 subtract ambient circuit 401. In other hardware configurations, the ambient light is removed 
from the red and infrared signals prior to the analog-to-digital converter. In those 
configiu^ations, the preprocessed input data can be directly sent to process red and infrared 
component circuit 203. The resultant signal is then transmitted to process red and infixed 
component circuit 203. The measured ambient components are also transmitted to ambient 

20 analysis circuit 402 where the characteristics of the ambient signal are used to help determine 
whether there is a probe off condition or excessive noise in the input signal. This ambient 
analysis is forwarded to data quality analysis circuit 208 which performs the error determination 
function for the pulse oximeter circuit to identify any aberrant conditions in the input signals 
which question the validity of the measurements taken as described below. , 

25 The red and infrared pulsatile components with ambient removed are applied to process 

red and infrared component circuit 203. The processing of the red and infixed components 
consists of several possible implementations illustrated in Figure 5. One implementation 
adaptively differentiates 501 the red and infrared signals to perform a digital differentiation of 
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these two signals to de-exponentiate the pulsatile component. Adaptive averaging circuit 502 
performs a local average of the red and infrared signals to estimate the non-pulsatile component. 

Alternatively, the differentiation can be replaced by logarithmic difference circuit 503 
which performs a digital subtraction of the logarithms of two sequential red signals and the 
5 same for two sequential infrared signals. The logarithmic differences are used in place of the 
differential arid local average values computed by circuits 501 and 502 respectively. 

As a result of these computations, a plurality of sets of processed red and infrared data 
values are stored in memory of data processing circuit 107 wherein each set consists of the 
processed red and infrared data values computed at two successive sampling intervals. 
10 Successive sets of data values therefore represent successive processed samples of 
measurements taken in successive time intervals. 

Whichever ones of these processes are selected, the ou4>ut processed signals are applied 
to compute "pulsatility index" circuit 204 which computes a "pulsatility index" indicative of the 
perfusion of the hemoglobin using a weighted combination of the red and infrared processed 
15 values. The perfusion is a term used to indicate the volume of arterial blood flow for a given 
volume of tissue and by looking at the red and infrared signals separately, an estimation of the 
present blood volume change in the tissue can be determined. The ""pulsatility index^' is 
therefore only a figure that can be correlated to the true perfusion of the tissue. 

These same processed signals output by process component circuit 203 are applied to 
20 regression computation circuit 205 which is the first step used in calculating the partial 
oxygenation Sp02 of the hemoglobin. The regression computation circuit 205 consists at least 
one of the plurahty of subcircuits illustrated in Figures 6 and 7. These circuits illustrate various 
combinations of regression coniputation and normalization. Circuit 601 is typically used to 
normalize the differentials received from process component circuit 203. The normalization 
25 circuit 601 multiplies each red differential by the corresponding inlfrared non-pulsatile intensity 
and multiplies each infrared differential by the corresponding red non-pulsatile intensity. The 
normalized differentials are then applied to regression circuit 602 which performs a linear 
regression fit using the normalized infrared differential as the independent variable and the 
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normalized red difTerential as the dependent variable. The slope of the best fit line (Figure 9) is 
a constant value that represents the best fit RRatio for that data set. A robust regression process 
can be implemented by removing the outliers in the collection of data points to thereby improve 
the least squares fit, or by other standard methods for robust linear regression. 

5 An alternative method of processing (shown in Figure 7) performs a regression fit to the 

differential values that are input to the regression computation circuit 205 and the output of the 
regression circuit is then normalized by normalize Slope circuit 603. In either case, the output 
of regression computation circuit 205 is the RRatio for the data set as well as a set of 
computation products that are forwarded to compute correlation coefficient circuit 206. In 
10 addition, the linear regression can be either a first order polynomial fit or a higher-order 
polynomial fit. In general, good input data should produce a fitted line having a finite slope, 
zero intercept and all higher order terms should be essentially zero. 

The compute correlation circuit 206 uses the RRatio calculated by regression 
computation circuit 205 as well as a set of parameters used in the computation of the RRatio to 

15 determine the goodness of the linear fit of this straight line in the set of red and infrared 
differentials that were transmitted by process component circuit 203 to regression computation 
circuit 205. A correlation coefficient is calculated which represents the goodness of the Unear 
fit of the straight line to the set of obtained data values. The intercept of the best fit line can be 
also computed to evaluate data quality. The correlation coefficient and zero offset values are 

20 transmitted to data quality analysis circuit 208 which uses these values as well as the ambient 
analyses obtained from process ambient input data circuit 202 to verify the quality of the input 
data. 

Probe-OfT Detection 

A probe-off condition can be detected by measuring tiie goodness of the linear fit* This 
25 means that the product of the slope and its inverse should equal one if the linear fit is a good 
one. 
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There is a standard statistical technique for measuring the goodness of a linear regression 
iit. The Fisher Linear Correlation essentially measures the goodness of the fit by using the 
constraint that the slope and its inverse must be a product very close to ±1 if the fit is good. . If 
the product of the slope and its inverse are significantly different than ±1, the data is bad. A 
5 higher order polynomial fit can be used instead of a linear and then the size of the higher order 
coefficients are used as a test. 

Data quality analysis circuit 208 monitors the correlation coefficient as well as other 
indicators to identify excessive input signal noise. A neural network or knowledge-based 
system can be used to identify a subset of anomalies indicative of errors. These anomalies can 

10 include n successive samples that exceed a predetermined variance and/or excessive variance in 
the ambient signal. Furthermore, the Sp02 varies slowly as a result of physiological changes 
and excessive correlation of changes in the determined Sp02 with changes in the ambient is an 
indication of an error condition. This circuit may also receive input pulse signal data available 
fit>m electrocardiograph machines which could be used to check the validity of the data and that 

1 5 the pulsatile rate of the data can be compared with the ECG values to identify motion artifacts in 
the input data. By the use of pattern matching, a feature vector can be computed to identify not 
only the presence of an error condition. It may also be used to classify the error. 

The best fit RRatio output by regression computation circuit 205 is transmitted to 
compute Sp02 circuit 207 where the partial oxygenation of the hemoglobin in the arterial blood 
20 is computed using an empirical formula to identify the correspondence between the input data 
and the actual oxygenation value as noted above. The results of this computation are then 
applied to display driver 109 for numeric output on display 115. While a specific embodiment 
of this invention has been disclosed, it is expected that those skilled in the art can and will 
design alternate embodiments of this invention that fall within the scope of the appended claims. 
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